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C§(Q 2 ) axial form factor is extracted from the ANL neutrino-deuteron 
scattering data with deuteron structure effects taken into consideration. 
The best fit of the C^(Q 2 ) axial form factor is obtained assuming dipole 
parametrization with C^(0) = 1.13 ± 0.15 and M A = 0.94 ± 0.08 GeV. 

PACS numbers: 25.30.Pt, 13.15. +g 

1. Introduction 

A knowledge of the cross sections for single pion production in neutrino- 
nucleon scattering is an important ingredient of the long base-line oscillation 
experiment analyses (see experiments: K2K [1] and T2K [2]). In particu- 
lar, prediction of the cross sections for l7r production in neutral current 
neutrino-nucleon scattering has great importance for the estimation of the 
background for observation of the — > v e neutrino oscillation. The reason 
is that the tt° decay signal into two photons can be mislead in the detector 
with the electron shower. 

More than twenty years ago two bubble chamber experiments: 12ft de- 
tector at Argonne National Laboratory (ANL) [3J 0] and 7ft detector at 
Brookhaven National Laboratory (BNL) [5] had collected enough data to 
investigate the W, Q 2 and energy dependence of the differential and total 
cross sections for neutrino-deuteron scattering. The ANL and BNL experi- 
mental measurements can still serve as an important source of information 
about cross sections for quasi-elastic and single pion production in neutrino 
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scattering. The nuclear effects (deuteron structure corrections) are relatively 
easy to consider, therefore analysis of these data allows to study neutrino- 
nucleon reactions. 

In this paper a part of the re-analysis [6] of the neutrino-deuteron scatter- 
ing data is presented. We re-examine the single pion production in neutrino 
deuteron interactions. Our aim is re-extraction of the C§(Q 2 ) axial form 
factor from the Itt data. In this article we discuss only the data collected 
at Argonne National Laboratory. In the longer contribution simultaneous 
analysis of both ANL and BNL data is described. 

The presentation is organized as follows in Sect. [2]the model for A(1232) 
excitation is shortly introduced. Sect. [3] contains description of the ANL 
data and details of the statistical analysis. Sect. |4] includes our final results 
and conclusions. 



2. Single pion production 

We consider the following process: 

u{k) + p{p) -> yT{k') + A++(p'){^ 7T+(0 +p(r)}. (1) 

By the k and k' the neutrino and muon four-momenta are denoted. The 
incoming nucleon momentum is denoted by p, while I and r are the final 
pion and proton momenta. The four momentum transfer is defined as follows 
Q 2 = —{k — k') 2 . The total hadronic momentum is p' = I + r. 

One of the simplest way to describe above reaction is to apply the 
Adler-Rarita-Schwinger formalism [7j. In this description the scattering 
amplitude is given by the contraction of the lepton with the hadronic 

(A ++ ,p'\ J^ c \p) currents. The hadronic current has Vector-Axial struc- 
ture, where the vector contribution is modeled by three unknown vector 
form factors C\ , C\ and C§. For these form factors we apply recent fits 
from Ref. 0. 

In general the axial contribution to the hadronic current is given by four 
form factors, however, additional constrains, motivated by Adler model |9j 
and PCAC hypothesis, reduce the number of unknown functions to the one: 
C£(Q 2 ). 

In the present analysis we consider the dipole parametrization of the 
C^(Q 2 ) axial form factor: 

Ci{Q 2 ) = C^)[l + ^ . (2) 
The axial mass Ma is a parameter to fit. We distinguish two fit cases: 
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(i) the value of C^(0) is established by PC AC and equals 1.15 [TO] : 

(ii) 0^(0) is treated as a free parameter. 

3. Re-analysis of the ANL data 

A subject of our analysis is the differential cross section data 

<I(J (Q?) = VezM) (3) 



dQ 2 

which was published in Ref. [4j. In the reaction (jTJ) the dominant contribu- 
tion comes from the excitation of the A(1232) resonance. The nonresonant 
background can be neglected. 

The proper analysis of the JTJ process requires taking into consideration 
the deuteron structure effects. In this work we follow the approach proposed 
in Ref. [IT]. From the practical point of view, in order to obtain the differen- 
tial cross section for neutrino-deuteron scattering, we multiply neutrino-/ree 
nucleon cr(Q 2 ) formula by 

R{Q 2 ) = <j{vd — » fi~ir + p + n)ja{yp — > fi^ir + p) 



correcting function. It is obtained from Ref. [TT] (for more details see. Ref. 
®)- 

To analyze the data we apply the standard x 2 approach. Besides statis- 
tical and non-correlated systematical uncertainties in each Q 2 bin, we take 
into account the overall uncertainty of the neutrino beam. It introduces to 
the standard x 2 formula additional normalization term: 




x 2 = V| -"•"•„,: """ ) --{'—r 1 ) ■ (4) 

where (J t h{Q 2 ) and a ex are the theoretical and experimental values of the 
differential cross sections. Acij denotes error in each Q 2 bin. r is the un- 
certainty of the neutrino flux and p is the normalization parameter which is 
going to be fitted. 

In order to compute the theoretical values of the cross section the follow- 
ing kinematical cuts are imposed: for the neutrino energy E £ (0.5, 6) GeV, 
and for the hadronic invariant mass W < 1.4 GeV. The data covers the 
range in Q 2 from Q 2 = 0.01 GeV 2 to Q 2 = 1 GeV 2 . 

The theoretical formula for the differential cross section in a given Q 2 
bin is the following: 



Ku P dQ 2 [ E up dE . . f 1A GeV 



°th(Q' Z ) = I ;'~"^2 / ' ^( E ) I dWa th (E,Q 2 ,W) (5) 
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M A [GeV] 


C 5 A (0) 


P 


X 2 /NDF 


free target 


0.948 ± 0.074 




1.15 ± 0.09 


1.50/7 


free target 


0.939 ± 0.082 


1.04 ± 0.14 


1.02 ± 0.19 


0.94/6 


deuteron 


0.937 ± 0.075 




1.03 ± 0.09 


0.81/7 


deuteron 


0.936 ± 0.077 


1.13±0.15 


1.02 ± 0.19 


0.80/6 



Table 1. The obtained numbers for C^(Q 2 ) and Ma resulting from the fitting 
procedure. 



with Q\ l0 = Q\- AQ?/2, Q\ uv = Qf + AQf/2. 

= / dE<f> ANL (E), (6) 

where E[ Q = 0.5 GeV and E up = 6 GeV. 

As it was discussed in Ref. [4] the normalization uncertainty of the total 
cross section, due to the lack of knowledge about the flux is estimated to 
be 15% for E € (0.5, 1.5) GeV and 25% for E > 1.5 GeV. Therefore in our 
discussion, we assume the average overall normalization uncertainty to be 
20%, r = 0.20. 

4. Results and conclusions 

In Tab. [T]the results of our fitting procedure are shown. In the analysis 
we consider four different cases. We start the fitting procedure without 
accounting deuteron correction. The results are shown in the first two rows 
of the Tab. [T]). Then the deuteron corrections are taken into consideration 
(look at the third and fourth rows of the Tab. []]). 

The comparison of our fits with the ANL data is shown in Fig. [Q where 
both the differential and total cross sections are plotted. 

Our final result (with deuteron effects) is : 



M A = 0.936 ± 0.077 GeV, C^(0) = 1.13 ± 0.15. 



(7) 
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Fig. 1. The differential (left panel) and total (right panel) cross sections for v+d — > 
fi~ +p + n + A ++ (1232) reaction measured at ANL experiment. The experimental 
points are taken from [4], The solid line denotes cross section computed with 
C5 (0) = 1.13 and Ma = 0.936 GeV. The dashed line represents cross section 
computed with C£(0) = 1.15 (fixed) and M A = 0.937 GeV. Both fits were obtained 
with accounting deuteron structure effects. The cut W < 1.4 GeV on hadronic 
invariant mass is imposed. The normalization parameter is p = 1. In the left panel 
in order to compute cross section the deuteron structure correction is imposed. 



As it can be noticed for all fits the obtained values of x 2 /NDF (NDF 
- number of degrees of freedom) are very good. We observe that the ax- 
ial mass, which is responsible for the shape of da/dQ 2 , is not affected by 
the deuteron corrections. The deuteron nuclear effects affect mainly the 
normalization of the cross section, therefore they are compensated by the 
normalization parameter p. 

The largest suppression of the a(Q 2 ), due to deuteron effects appears in 
low-Q 2 (below 0.1 GeV 2 ). That is why the fit of C^(0) without deuteron 
correction is relatively low (rs 1.04). When the deuteron effects are ac- 
counted the obtained value is « 1.13 which is very close to the PCAC value. 
However, both results for C§(0) are always in agreement with the PCAC 
value. 

The fits obtained for the ANL data turns out to be quite similar to those 
obtained from the simultaneous fit to both the ANL and BNL data. The 
full presentation of this analysis will be published in Ref. [6]. 
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